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ABSTRACT 

An  analytical  study  was  undertaken  to  determine  the 
two-dimensional  wall  conduction  effects  in  an  internally 
finned,  rotating  heat  pipe.  The  finite  element  method  was 
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with  triangular  fins.  Heat  transfer  rates  were  shown  to  be 
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I .  INTRODUCTION 


A.  THE  ROTATING,  NON -CAP ILLARY  HEAT  PIPE 

The  rotating,  non-capillary  heat  pipe  is  a  self-contained 
device  designed  to  transfer  large  quantities  of  heat  in  a 
rotating  system.  The  heat  pipe  contains  an  evaporator,  a 
condenser,  and  a  working  fluid  as  shown  in  Figure  1.  The 
working  fluid  is  vaporized  by  heat  addition  at  the  evaporator. 
The  vapor  flows  toward  the  condenser  end  as  a  result  of  a 
vapor  pressure  differential.  There,  the  vapor  condenses  on 
the  inner  wall,  releasing  its  latent  heat  of  vaporization. 

The  condensate  then  returns  to  the  evaporator  where  the  cycle 
is  repeated. 

The  difference  between  the  conventional  heat  pipe  and  the 
rotating,  non-capillary  heat  pipe  is  the  condensate  return 
mechanism.  A  conventional  heat  pipe  employs  the  capillary 
pumping  action  of  an  internal  porous  wick  to  transport  the 
condensate  to  the  evaporator.  The  disadvantage  of  this  scheme 
is  the  resistance  to  condensate  flow  through  the  wick.  Addi¬ 
tionally,  the  presence  of  nucleate  boiling  in  the  evaporator 
could  create  bubbles  inside  the  wick  structure,  reducing  the 
effectiveness  of  the  pumping  action  [1]^.  In  the  rotating 
heat  pipe,  rotation  about  its  longitudinal  axis  establishes 
a  centrifugal  force  field  which,  when  combined  with  a  slight 

^  Numbers  in  brackets  indicate  references  listed  in  the 
Bibliography . 
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taper  of  the  condenser  section,  forces  the  condensate  toward 
the  evaporator.  At  high  rotation  rates,  the  heat  transfer 
rate  is  significantly  increased  over  the  rate  achieved  in 
conventional  heat  pipes  of  similar  dimensions. 

B .  BACKGROUND 

Previous  analyses  [2-9]  of  the  rotating  heat  pipe  have 
been  conducted  to  determine  its  performance  characteristics. 
These  studies  determined,  among  other  things,  that  under 
normal  circumstances,  the  heat  transfer  rate  is  limited  pri¬ 
marily  by  the  amount  of  heat  which  can  be  transferred  through 
the  condenser  wall. 

In  an  effort  to  improve  the  condenser  heat  transfer  per¬ 
formance,  Schafer  [8]  estimated  and  compared  the  thermal 
resistances  across  the  condensate  film,  the  wall,  and  the 
outside  convective  film  using  Woodard's  earlier  data  [9]. 

He  showed  that  with  the  use  of  outside  spray  cooling  and  a 
thin  copper  wall,  the  controlling  resistance  is  the  inside 
film  resistance.  Schafer  suggested  that  this  resistance  can 
be  reduced  (1)  by  the  proper  selection  of  working  fluid,  and 
(2)  through  the  use  of  internal  axial  fins  as  sketched  in 
Figure  2a. 

Machining  axial  fins  into  the  inner  surface  of  the  con¬ 
denser  would  improve  heat  transfer  performance  through  two 
mechanisms:  a  decrease  in  condensate  film  thickness  near 

the  apex  of  each  fin,  and  an  increase  in  inner  wall  surface 
area.  However,  with  the  addition  of  fins,  the  wall  thermal 
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resistance  would  also  increase  due  to  the  increase  in 
effective  wall  thickness.  Of  course,  the  use  of  copper  as 
the  wall  material,  because  of  its  high  thermal  conductivity, 
would  minimize  the  added  wall  resistance. 

Schafer  then  developed  an  analytical  model  of  the  intern¬ 
ally  finned  condenser  for  the  case  of  fins  of  triangular 
profile  (Figures  2a,  2b).  In  the  analysis,  he  assumed  one¬ 
dimensional  heat  conduction  through  the  wall  and  assumed  a 
linearly  decreasing  temperature  profile  from  the  fin  apex 
(point  A)  to  the  midpoint  of  the  trough  (point  C) .  Results 
from  the  analysis  indicated  significantly  improved  heat 
transfer  rates  for  the  finned  condenser  compared  to  the  smooth 
condenser.  He  pointed  out,  however,  that  two-dimensional 
conduction  effects  may  be  significant  so  that  a  more  detailed 
analysis  of  the  heat  conduction  problem  within  the  wall  was 
necessary . 

C.  OBJECTIVES 

The  objectives  of  this  study  were,  therefore,  to  (1) 
develop  an  analytical  solution  to  account  for  two-dimensional 
conduction  effects  within  the  wall  for  the  triangular  fin 
profile,  and  (2)  compare  this  solution  with  the  solution 
derived  by  Schafer. 
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II.  THEORETICAL  PROGRAM 


A.  REVIEW  OF  SCHAFER  ANALYSIS  [8] 

Schafer's  study  considered  an  internal  fin  geometry 
represented  by  Figure  2a.  Since  the  sum  of  wall  thickness 
and  fin  height  was  assumed  to  be  much  smaller  than  condenser 
radius,  an  average  radius  from  any  point  on  the  axis  to  the 
condenser  wall  was  taken  to  be  the  distance  from  the  axis  to 
a  position  midway  between  the  base  of  the  fin  and  the  apex 
of  the  fin,  r^ .  The  centrifugal  acceleration  term  was  then 
evaluated  at  this  average  radius,  r^w^. 

The  thin  wall  assumption  also  allowed  the  curvature  of 
the  wall  to  be  neglected.  The  problem  could  then  be  analyzed 
using  a  planar  conduction  heat  transfer  model  as  shown  in 
Figure  2b.  Furthermore,  because  of  repetitive  geometry  the 
problem  could  be  reduced  to  the  simpler  geometry  shown  in 
Figure  3.  Cartesian  coordinates  were  then  selected,  as  shown, 
to  permit  a  simplified  fin  condensation  analysis. 

Assumptions  in  the  analysis  were  based  upon  the  earlier 
work  of  Ballback  [2] .  These  assumptions  were  similar  to  those 
used  in  the  Nusselt  analysis  of  film  condensation  on  a  wall. 
The  more  important  of  these  assumptions  are; 

(1)  steady  state  operation, 

(2)  film  condensation,  as  opposed  to  dropwise  condensation, 

(3)  laminar  condensate  flow  along  both  the  fin  and  the 
trough, 

(4)  static  balance  of  forces  within  the  condensate. 
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(5)  one- dimensional  conduction  heat  transfer  through  the 
film  thickness  (no  convective  heat  transfer) , 

(6)  no  liquid-vapor  interfacial  shear  stresses,  and 

(7)  no  condensate  subcooling. 

Additional  assumptions  in  the  analysis  are  noted  later. 

The  complete  development  of  Schafer's  finned  condenser 
heat  transfer  theory  is  presented  in  his  thesis  [8] .  However, 
since  the  fluid  flow  and  heat  transfer  theory  was  applicable 
to  the  two-dimensional  conduction  heat  transfer  analysis,  a 
summary  of  his  work  is  presented  here,  prior  to  the  develop¬ 
ment  of  the  two-dimensional  conduction  model. 

1 .  Fin  Condensate  Mass  Flow 

A  Nusselt  analysis  was  performed  for  film  condensation 
on  the  wall  of  the  fin.  In  this  analysis,  the  fluid  momentum 
equation  was  constructed  to  account  for  the  effects  of  condenser 
cone  half- angle,  (p ,  fin  half -angle,  a,  and  the  centrifugal 
force  field.  The  pressure  differential  term  was  neglected. 
Condensate  flow  on  the  fin  was  assumed  to  be  in  the  z-direction 
only. 

An  expression  for  fluid  velocity,  w(y,z),  was  then 
derived  by  substituting  the  shear  stress  equation  for  laminar 
flow,  ^  ,  into  the  momentum  equation  and  solving 

for  w(y,z).  Average  velocity,  w(z),  was  then  found  by  inte¬ 
grating  w(y,z)  over  y  from  y  =  0  to  y=5(z),  the  local  film 
thickness.  This  average  velocity  was  then  placed  in  the  fin 
mass  flow  equation  to  yield  an  expression  for  fluid  mass  flow 
rate  along  the  fin  per  unit  condenser  length  as  a  function 
of  6 (z)  , 
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(1) 


Mj(z)  = 


P£  co^CR^  +  X  sin  ({))  cos  (i)  cos  a  <S^(z) 


3y. 


2 .  Trough  Condensate  Mass  Flow  Equation 

Fluid  flow  within  the  trough  was  assumed  to  be  one¬ 
dimensional  in  the  x-direction  and  of  constant  velocity  across 
the  width  of  the  trough.  Therefore,  condensate  flowing  from 
the  fins  into  the  trough  was  assumed  to  accelerate  immediately 
to  the  velocity  of  the  trough  fluid.  Cross-sectional  area  of 
the  trough  condensate  region  was  twice  the  area  shown  in 
Figure  4. 

The  x-direction  momentum  equation  for  an  infinitesimal 
fluid  element  in  the  trough  was  constructed  with  the  pressure 
differential  term  neglected.  The  trough  mass  flow  expression 
was  then  developed  in  a  manner  similar  to  the  development  of 
the  fin  mass  flow  equation  (1) .  The  average  trough  velocity 
was  derived  and  substituted  into  the  trough  mass  flow  equation 
to  yield  total  mass  flow  through  the  trough  region, 


+x  sin  sin  <(> 


(£6*Cx)  +  sin  a  cosar5*(x)). 


C2) 


where  the  expression  in  the  second  parenthesis  is  the  cross- 
sectional  area  of  the  trough  condensate  region. 

Total  condensate  mass  flow  was  identical  to  trough 
condensate  mass  flow  since  all  condensate  eventually  flowed 
to  the  trough.  Therefore,  once  the  total  mass  flow  rate  at 
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any  position  x  was  known,  trough  condensate  thickness,  6*(x), 
could  be  determined  using  equation  (2). 

3 .  Condenser  Heat  Transfer  Determination 

Schafer  assumed  zero  heat  flux  boundary  conditions  on 
both  sides  of  the  wall  section  shown  in  Figure  4.  Temperature 
at  the  fin  apex  was  taken  to  be  fluid  saturation  temperature 
since  film  thickness  at  this  position  was  assumed  to  be  zero. 

He  then  imposed  a  linearly  decreasing  temperature  on 
the  inner  wall  surface  from  the  apex  of  the  fin  to  the  mid¬ 
channel  position.  The  surface  temperature  expression  became 

Ts  -  T^(z)  =  32,  C3) 

where  3  was  an  unknown  constant.  This  choice  of  temperature 
distribution  was  arbitrary  but  allowed  a  simplified  one¬ 
dimensional  conduction  analysis  of  the  condenser  wall.  He 
then  imposed  the  restriction  that  no  heat  transfer  would 
occur  through  section  BE  of  the  wall  (Figure  3) .  All  heat 
transfer  through  the  fin  surface,  thus,  flowed  out  of  the 
wall  through  section  FE. 

Heat  transfer  through  the  fin  was  equated  with  one¬ 
dimensional  heat  transfer  through  a  rectangular  geometry  of 
the  same  width  as  section  FE,  and  of  height,  t*,  equal  to  t 
plus  one-half  the  fin  height  (Figure  4) .  Temperature  of  the 
upper  surface  of  the  rectangle  was  set  equal  to  the  tempera¬ 
ture  at  a  point  midway  between  positions  A  and  B.  The 
expression  for  heat  transfer  rate  per  unit  condenser  length 
through  the  wall  and  outer  convective  thermal  resistance 
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of  the  fin  became 


t  CTj  -  )  -  T„  ] 


•  * 

Q^(z3 


t* 


‘out 


(4) 


It  was  then  desired  to  develop  an  expression  for 

•  * 

Qj  (z)  in  terms  of  the  heat  transfer  rate  through  the  fin 
condensate  film.  First,  the  equation  for  fin  condensate 
film  thickness  was  derived  by  differentiating  (1)  with  respect 
to  z  and  substituting  the  result  into  the  differentiated 
equation  for  mass  flow  rate  per  unit  condenser  length. 


The  expression  produced  was  algebraically  rearranged  to  form 
a  differential  equation  in  6(z), 


6^(z)d6(z)  = 


kf  (Ts  -  Pf  dz 


P£  0)^  (R^  +  X  sin  (J) )  h£g  cos  ({)  cos  a 


(6) 


Equation  (3)  was  then  substituted  in  (6)  and  the  result 
integrated  over  z  from  z  =  0  to  z  to  give  an  equation  for 
6(z), 


<Slz)  =  r- 


2  k£  6  P£ 


p/  (R^  +  X  sin  4>)h  £g  cos  4)  cos  a 


1  V 

z  (7) 


The  value  for  fin  heat  transfer  rate  per  unit  condenser 
length  was  then  evaluated  using 

Q^(z*)  =  h£g  M^(z*)  (8) 
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where  z*  is  defined  in  Figure  4.  This  expression  was  set 
equal  to  (4)  to  yield  a  transcendental  equation  in  3, 


[(T. 


3Zo 


-  T  ] 

QD  J 


^O 

T* 


out 


1/  V4  3/ 

(2K^g)V‘*  (p^oj)  '^  [  (R  +x.sin  <)))  h^  cos  (j)  cos  a]  z*  ^ 


(9) 


where  was  the  width  of  the  base  of  the  fin.  The  expressions 
on  the  left  and  right  hand  sides  of  this  equation  were  equiva¬ 
lent  only  when  all  boundary  conditions  were  satisfied.  The 
task  was  then  to  solve  iteratively  for  3  until  equation  (9) 
was  satisfied.  (z*)  was  then  determined  in  equation  (4) 
using  the  final  iterated  value  of  3. 

Heat  transfer  rate  from  the  trough  region  was  evaluated 


using 


Q^Cx) 


(T  -  T  ) 
^  s  avg-^ 

6*(x) 


(10) 


where  calculated  as  an  average  temperature  of  the 

trough  wall  and  e  was  the  total  trough  width  (Figure  4) . 

A  flowchart  description  of  the  computer  program  used 
in  the  Schafer  analysis  is  given  in  Figure  5.  The  axial 
length  of  the  condenser  was  first  divided  into  increments 
of  equal  length,  Ax.  Solution  of  the  heat  transfer  problem 
was  then  started  by  assuming  an  initial  trough  condensate 
thickness,  6*^  =  0.00006752  feet.  This  initial  value  was 
based  upon  a  mass  flow  analysis  by  Sparrow  and  Gregg  of 
condensation  on  a  rotating  disk  [10]  .  It  was  assumed  in 
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Schafer's  analysis  that  the  effect  on  the  Sparrow  and  Gregg 
solution  of  enclosing  the  disk  within  the  condenser  would 
be  minimal.  A  study  of  the  effects  of  a  ten- fold  increase 
or  decrease  from  the  above  value  of  6*^  indicated  that  the 
initial  6*^  assumption  would  not  affect  the  heat  transfer 
solution  by  more  than  ten  percent.  Therefore,  the  deter¬ 
mination  of  an  exact  value  of  initial  trough  condensate 
thickness  was  not  crucial  to  the  problem. 

Calculation  of  condensate  mass  flow  contribution  from 

the  endplate  then  followed.  Then,  for  each  increment  Ax,  an 

3  z 

initial  value  of  3  was  calculated  by  removing  the  term 

from  the  left  side  of  equation  (9)  and  solving  directly  for 
3.  The  full  equation  was  then  employed  to  solve  iteratively 
for  3.  Calculation  of  incremental  fin  heat  transfer  rate 
then  followed  by  multiplying  the  value  obtained  in  (4)  by 
the  incremental  length.  Ax.  Incremental  trough  heat  transfer 
rate  was  next  evaluated  by  multiplying  equation  (lO)  by  Ax. 
Incremental  trough  and  fin  heat  transfer  rates  were  then 
summed  to  yield  the  total  incremental  heat  transfer  rate  for 
the  geometry  of  Figure  4. 

Total  incremental  mass  flow  rate  was  evaluated  by 
dividing  twice  the  sum  of  the  incremental  trough  and  fin  heat 
transfer  rates  by  the  latent  heat  of  vaporization  of  the  fluid. 
As  the  incrementation  progressed,  incremental  mass  flow  rates 
were  added  to  mass  flow  from  preceding  increments  to  yield 
C^) •  Equation  (2)  was  then  solved  for  6*(x)  with  the 
aid  of  a  polynomial  rootfinder  subroutine.  The  new  value  of 
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6*(x)  was  then  used  to  begin  the  solution  for  the  next 
increment.  When  the  total  length  of  the  condenser  had  been 
traversed,  total  heat  transfer  rate  was  evaluated. 


^tot 


=  2  N 


fin 


100 

I 

i=l 


Ax[Q  (z*).  + 


Q,(x)i] 


B.  TWO-DIMENSIONAL  HEAT  TRANSFER  SOLUTION 


(11) 


1 .  Introduction 

In  this  study,  two-dimensional  conduction  through  the 
condenser  wall  was  considered  in  the  heat  transfer  analysis 
using  a  finite  element  solution.  The  finite  element  method 
was  chosen  over  the  method  of  finite  differences  because  of 
considerations  involving  wall  geometry  and  computer  utiliza¬ 
tion  time.  Myers  [11]  has  shown  that  when  these  methods  are 
compared  for  the  case  of  two-dimensional  conduction  in  a 
rectangular  block  with  heat  generation,  the  finite  element 
method  more  closely  approximates  the  exact  solution  when  the 
number  of  nodal  points  considered  in  each  method  are  equal. 
This  is  particularly  true  when  the  nodal  point  meshes  are 
coarse . 


To  achieve  the  same  accuracy  as  with  the  finite  element 
method,  the  finite  difference  method  must  use  a  finer  nodal 
mesh,  especially  for  non-rectangular  geometries.  Additionally, 
when  fins  of  different  shape  are  considered,  the  mesh  of  the 
finite  difference  model  must  be  completely  reconstructed  for 
each  new  shape.  In  the  finite  element  program,  however,  a 
simple  change  in  nodal  coordinates  is  all  that  is  necessary 
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to  accommodate  the  change  in  fin  shape.  Since  future  studies 
may  analyze  the  use  of  varied  fin  shapes  with  the  rotating 
heat  pipe,  it  was  desired  to  incorporate  the  more  flexible 
finite  element  method  in  the  computer  program  to  avoid  having 
to  rearrange  the  program  for  each  fin  shape. 

While  the  accuracy  of  the  finite  element  method  is 
superior,  the  finite  difference  method  requires  less  computer 
utilization  time  for  the  same  number  of  nodal  points.  If  a 
finite  difference  model  of  finer  mesh  was  chosen,  the  decrease 
in  computer  time  due  to  the  method  might  be  offset  by  an  in¬ 
crease  in  time  required  for  the  increased  number  of  nodes. 

The  overall  improvement  in  computer  utilization  time  might, 
therefore,  be  minimal.  As  a  result,  in  this  study,  the  choice 
of  method  was  not  based  upon  computer  time  requirements.  Also, 
a  finite  element  computer  program  [12]  was  available  which 
could  be  adapted  for  the  purposes  of  this  analysis.  For  the 
reasons  outlined  above,  therefore,  the  finite  element  method 
was  selected. 

The  steady  state  differential  equation  for  conduction 
heat  transfer  is  given  in  Arpaci's  text  [13]  as 


3  ,  ,  3T  s  ^  3  .  ,  3T  .  ^ 
^  (  k  )  +  ^rk  -rrrJ  O* 
3y  y  3y  3z  ^  t  a? 


"z  3  z 


(12) 


Boundary  conditions  are  normally  of  two  kinds;  those  of 
specified  temperature  on  a  boundary. 


T  =  T^  ,  (boundary  b  only) 


(13) 


and  those  of  specified  heat  flux  or  convective  condition 
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on  the  boundary, 


+  [h(T  -  T^)  ]  cijj  =  0 .  (boundary  b  only) 


(14) 


where  H  and  i  are  direction  cosines  o£  the  outward 
y  z 

normal  to  the  boundary's  surface. 


2.  The  Finite  Element  Method 


Many  references  including  [11],  [14],  and  [15]  develop 


the  concept  of  using  material  elements  of  finite  dimensions 
within  a  structure  to  model  the  differential  equation  and  its 
boundary  conditions.  The  steady-state  field  equation, (12)  , 
may  be  discretized  by  considering  that  the  temperature  at  any 
point  within  an  element  is  a  function  of  the  temperatures  of 
specified  nodes  on  that  element. 


T  =  <N>  {  T  }®, 


(15) 


where  <(N^  is  the  row  matrix  shape  function  of  position  and 
{ T  is  the  nodal  temperature  vector  for  the  element.  The 


field  equation  and  boundary  equations  for  each  element  may 
be  mathematically  rearranged  and  integrated  to  yield  a  matrix 
system  of  equations  for  that  element. 


[H]®  {  T  =  -{  , 


(16) 


where 


and 
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(18J 


1 

*  (/^JN)h<N>dc)^  -  MT„)^c^. 

Equation  (18)  is  applied  only  over  relevant  boundaries.  The 
element  matrix  equations  are  then  assembled  into  a  master 
matrix  equation  for  the  entire  structure.  Solution  by 
numerical  methods  then  follows. 

The  convective  boundary  conditions  were  applied  in 
the  finite  element  solution  of  the  finned  condenser  using  a 
lumped  convective  condition  formulation.  Average  nodal  con¬ 
vective  coefficients,  h  ,  were  determined  for  each  node 

avg 

subject  to  a  convective  environment.  The  nodal  h  values 

avg 

were  then  considered  to  be  constant  over  the  boundary  region 
associated  with  the  node.  The  method  of  determination  of 
these  coefficients  is  described  in  a  later  section. 

Equation  (18)  was,  thus,  reduced  to  the  following 

form , 

{F}  ®  =  ({h^,^„  T}  ®  T  C19) 

b  avg  b  avg  <»  b  b 

This  lumping  procedure  reduced  the  complexity  of  the  convective 
condition  calculations  and  did  not  seriously  affect  the  values 
within  the  master  matrix,  [H] . 

Reference  [15]  presents  a  development  of  the 
isoparametric  finite  element  concept.  This  concept  allows 
the  solution  of  equations  (12)  ,  (13)  ,  and  (14)  for  elements 
of  non-uniform  shape.  The  application  of  the  method  of  finite 
elements  is,  therefore,  extended  to  a  wider  range  of  physical 
problems . 
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The  eight-nodal  finite  element  mesh  used  in  the 
two-dimensional,  analysis  of  the  axially  finned  condenser  is 
shown  in  Figure  6.  A  computer  program  for  the  case  of  two- 
dimensional  conduction  was  developed  using  a  modified  version 
of  the  three-dimensional,  isoparametric  finite  element  program 
assembled  by  Lew  [12] . 

3 .  Test  Cases 

Before  the  finite  element  solution  was  applied  to  the 
condenser  problem,  a  series  of  test  cases  were  analyzed. 

These  cases  allowed  a  comparison  of  analytical  solutions  to 
results  obtained  using  the  finite  element  program.  The  two- 
dimensional  heat  conduction  problems  studied  were 

(1)  a  rectangular  block  with  specified  surface  temperature, 

(2)  a  rectangular  block  in  a  uniformly  convective  environ¬ 
ment, 

(3)  a  rectangular  block  with  a  triangular  fin  in  a  uniformly 
convective  environment, 

(4)  a  rectangular  block  with  a  convex  parabolic  fin  in  a 
uniformly  convective  environment,  and 

(5)  a  rectangular  block  with  a  concave  parabolic  fin  in  a 
uniformly  convective  environment. 

Results  of  the  test  cases  are  presented  in  Figures  7-11  with 
values  obtained  using  the  finite  element  program  placed  in 
brackets.  For  the  first  two  cases,  the  heat  transfer  rates 
calculated  with  the  aid  of  the  finite  element  computer  program 
agreed  well  with  the  exact,  two-dimensional  solutions. 

Although  finite  element  heat  transfer  solutions  were  rela¬ 
tively  accurate,  nodal  temperature  values  were  slightly 
different  from  the  exact  values.  This  temperature  discrepancy 
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led  to  numerical  difficulties  when  the  finite  element  solution 
was  applied  to  the  stainless  steel  condenser.  These  diffi¬ 
culties  will  be  discussed  more  fully  in  Chapter  III. 

For  the  finned  block  cases,  analytical  expressions 
were  derived  for  heat  transfer  rates  through  extended  surfaces 
using  the  one-dimensional  fin  analysis  as  described  in  [16] . 
Heat  transfer  rates  calculated  by  the  finite  element  program 
for  this  uniformly  convective  environment  did  not  vary  sig¬ 
nificantly  from  the  one-dimensional  results.  Since  test  case 
results  from  the  finite  element  program  agreed  favorably  with 
the  derived  results,  the  program  was  assumed  to  be  operating 
properly  and  was  applied  to  the  condenser  heat  transfer 
problem. 

4 .  Condenser  Boundary  Condition  Application 

The  axially  finned  condenser  problem  was  analyzed 
by  considering  the  geometry  of  Figure  6.  Boundary  conditions 
for  the  problem  were  considered  to  be  the  same  as  those  in 
the  Schafer  analysis  with  the  exception  of  the  inner  surface 
temperature  distribution.  Instead  of  using  a  linearly  decreas¬ 
ing  temperature  profile  along  the  wall  from  the  apex  of  the 
fin  to  the  midpoint  of  the  channel,  a  parabolic  distribution 
was  assumed  and  used  to  calculate  average  nodal  convective 
coefficients . 

The  Nusselt  analysis  of  the  fin  condensation  problem 
yielded  a  differential  equation  in  6(z),  fin  condensate 
thickness,  which  is  identical  to  equation  (6)  used  in  Schafer's 
analysis . 
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Consistent  with  the  eight-nodal  finite  element 
concept,  a  parabolic  temperature  distribution  was  then 
assumed  to  exist  along  the  fin  surface  using  the  boundary 
condition,  T  =  at  z  =  0,  since  6(0)  =  0.  Therefore, 

T^(z)  =  ^  b.z  *  T3  ,  (20) 

or 


Ts  -  =-a,z^  -  b,z  ,  (21) 

where  a^  and  b^  were  unknown  constants. 

When  (21)  was  substituted  into  equation  (6)  and  the 
result  integrated  from  z  =  0  to  z ,  the  expression  for  6(z) 
became 


6(z) 


p£^  w^(R  +  X  sin  4))  h£g  cos  (})  cos  a 


(22) 


The  local  convective  heat  transfer  coefficient,  assuming 
one-dimensional  heat  conduction  through  the  thin  condensate 
film,  was  then 


h(z) 


fin 


^f 

^) 


=  [■ 


k  ^ 
^f 


P£^  (R  +  X  sin  (})  )  hg^  cos<|)  cos  a 


4y£  (- 


Z  3 


-  b 


V4 

]  , 


(23) 


and  was  applicable  from  z  =  0  to  z  =  z*  =  z^  -  6*/cos  a 
(Figure  6)  . 


Across  the  width  of  the  trough,  from  z  =  z^  to  the 
mid-channel  position. 


k£ 

trough  6*(x) 


(24) 
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6*(x)  was  considered  constant  across  the  trough  width.  Based 
upon  a  check  of  the  wall  temperatures  generated  from  Schafer's 
computer  program,  the  fluid  temperature  difference  across  the 
width  of  the  trough  was  predicted  to  be  small.  Fluid  thermal 
conductivity  values  were,  therefore,  assumed  to  remain 
relatively  constant.  As  a  result,  was  assumed  to 

be  constant  across  the  trough  width  and  was  calculated  using 
equation  (24)  and  a  representative  average  inner  wall  temper¬ 
ature  to  determine  k£. 

An  h(z)  temperature  distribution  was  then  needed 

between  z  =  z*  and  z  =  z  .  Since  the  thickness  of  condensate 

o 

film  increased  linearly  from  the  value  at  z*  to  the  value  at 
z^,  the  convective  coefficient  was  assumed  to  decrease  linearly 
from  h(z*)f.j^  to  h(x)„ough- 

Other  boundary  conditions  were  then  applied  in 
accordance  with  the  assumptions  of  Schafer,  namely,  zero  heat 
flux  through  the  sides  of  the  geometry,  T  =  at  the  fin 
apex,  and  constant  convective  coefficient,  ,  along  the 

outside  surface  of  the  wall.  Boundary  conditions  for  the 
problem  were,  thus,  completely  specified. 

In  his  expression  for  trough  condensate  cross-sectional 
area,  Schafer  assumed  that  the  surfaces  of  the  condensate 
flows  from  the  fin  and  the  trough  were  coincident  at  a  position 
directly  above  the  trough-fin  wall  juncture  (Figure  4) .  For 
small  fin  film  thicknesses  and  other  than  small  fin  angles, 
this  assumption  was  not  accurate.  In  this  study,  it  was 
assumed  that  fin  condensate  thickness  was  much  smaller  than 
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trough  condensate  thickness.  The  trough  condensate  region 
could,  therefore,  be  more  accurately  described  as  shown  in 
Figure  6.  The  modified  trough  mass  flow  equation ^became 

P£^  (R^  sin  <}>)6*(x)^  sin  <p 

"  T  (6*(x)e  +  6*(x)2tana)  . 

:)Uf 

(25) 

It  should  be  noted  that  this  equation  differs  from  equation 
(2)  only  in  the  expression  for  cross-sectional  area  of  the 
trough  condensate  region.  The  change  in  cross-sectional  areas 
involved  the  term,  1/cos^  a,  which  for  large  fin  angles,  could 
make  a  significant  difference  in  cross-sectional  areas  of  the 
two  equations. 

5 .  Condenser  Heat  Transfer  Solution 

Within  the  computer  program,  the  heat  transfer  rate 
for  the  entire  condenser  was  determined  by  dividing  its  axial 
length  into  one  hundred  equal  length  increments  of  length,  ^x, 
solving  for  the  incremental  heat  transfer  rate  using  an  iter¬ 
ative  procedure,  and  summing  those  incremental  rates  to  yield 
a  total  heat  transfer  rate.  A  flowchart  description  of  the 
computer  program  used  in  this  analysis  is  given  in  Figure  12. 

At  each  increment,  constants  a^  and  b^  were  determined 

using  a  Lagrangian  interpolation  fit  of  the  temperatures 

/• 

calculated  from  the  previous  iterate  at  nodes  1,  3,  and  6. 

For  the  first  increment,  a^  and  b^  were  initially  calculated 
from  an  initial  guess  of  the  temperatures  at  those  nodes. 
Convective  coefficients  at  several  locations  along  the  fin 
were  determined  using  equation  (23)  and  the  linearly  decreasing 
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h(z)  assumption  from  z  =  z*  to  z  =  z^.  Next,  average  nodal 
values,  hj^^g,  and  the  average  h  o£  the  left  hand 

boundary  region  associated  with  node  6,  (Figure  13),  were  then 
calculated  using  the  position  dependent  values  of  hCz)^^^  in 
a  Simpson's  Rule  approximation.  trough  then  calculated 

from  equation  (24)  using  the  6*  value  calculated  from  the 
previous  increment.  This  value  was  used  as  the  average 
convective  coefficient  for  nodes  9  and  12  and  the  right  hand 
boundary  region  of  node  6.  The  average  convective  coefficients 
calculated  for  the  left  and  right  hand  boundary  regions  of 
node  6  were  then  weighted  by  their  respective  region  lengths. 
The  two  values  were  then  summed  to  yield  the  average  convective 
coefficient  for  node  6,  hg^^g.  The  remaining  boundary  condi¬ 
tions  were  then  applied  and  incremental  heat  transfer  rate 
per  unit  condenser  length,  ,  and  nodal  temperature  values 
were  generated  from  the  finite  element  subroutines. 

If  the  resultant  heat  transfer  rate  converged  to 
within  the  criteria  of  0.05  percent  change  from  the  value  of 
the  previous  iterate. 


<  0.0005  , 


(26) 


the  incremental  problem  was  deemed  solved  and  set  equal  to 
the  final  iterate,  Qj .  If,  on  the  other  hand,  the  convergence 
criteria  was  not  met,  the  newly  calculated  nodal  temperature 
values  were  used  to  establish  a  new  set  of  inner  wall 
convective  coefficients.  After  the  remaining  boundary 


29 


conditions  had  been  applied,  the  finite  element  solution  was 
repeated  until  heat  transfer  rate  converged  to  within  the 
acceptable  0.05  percent  tolerance.  Because  Qj  was  initial¬ 
ized  at  zero  for  each  increment,  the  iteration  procedure  was 
repeated  at  least  twice  for  every  increment. 

Incremental  mass  flow  rate  was  then  calculated  using 


A  M 


tot. 


(27) 


This  value  was  then  added  to  the  mass  flow  rate  from  previous 
increments.  The  subsequent  interval's  trough  condensate 
thickness  was  then  calculated  using  a  polynomial  root-finder 
subroutine  to  solve  equation  (25) .  The  incrementation 
continued  until  the  entire  length  of  the  condenser  had  been 
traversed.  Incremental  heat  transfer  rates  were  than  summed 
and  a  total  heat  transfer  rate  calculated  from  (3) . 

100  • 

•^tot  “  2  Nf  (Q.  AX)  ,  (28) 

where  was  the  number  of  fins  along  the  inner  wall 

circumference . 

To  start  the  incrementation,  two  initial  assumptions 

were  made.  The  initial  increment's  value  of  5*  =  6*  was 

1 

taken  from  the  analysis  of  Schafer  as  0.00006752  feet.  For 
reasons  explained  in  the  previous  sub-chapter,  a  large  variance 
in  the  accuracy  of  this  value  would  not  significantly  affect 
the  heat  transfer  results  obtained  from  the  computer  program. 
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Additionally,  at  the  first  increment,  initial  guesses 
of  fin  temperatures  at  nodes  1,  3,  and  6  were  made  in  order 
to  begin  the  solution.  As  an  example,  the  temperatures  used 
to  begin  the  solution  for  the  copper  condenser  with  fluid 
saturation  temperature  of  200®F  were  =  200®F,  =  195°F, 

and  T  =  191°F. 

6 

For  the  copper  condenser,  the  use  of  widely  variant 
initial  nodal  temperature  values,  =  200°F,  =  191°F, 

and  T  =  182®F  had  little  effect  on  convergence  to  the  heat 
transfer  solution  for  the  first  increment.  However,  as  will 
be  discussed  in  Chapter  III,  the  initial  temperature  guess, 
coupled  with  the  slight  inaccuracy  of  the  finite  element 
nodal  temperature  solution  produced  numerical  difficulties 
within  the  computer  program  for  the  stainless  steel  condenser. 
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III.  DISCUSSION  OF  RESULTS 


Computer  results  for  the  finned  copper  condenser  were 
generated  for  three  outside  convective  coefficients,  three 
fin  angles,  and  five  rotation  rates.  The  working  fluid  used 
throughout  the  analysis  was  water.  Representative  results 
are  presented  in  Figures  14  and  16,  17,  18. 

Figure  14  compares  the  heat  transfer  results  of  this 
study  with  the  results  of  the  analysis  of  Schafer  for  given 
operating  conditions.  It  can  be  seen  that  the  two-dimensional 
conduction  analysis  allowing  a  parabolic  temperature  distri¬ 
bution  along  the  fin  predicted  heat  transfer  rates  for  the 
finned  condenser  which  are  approximately  seventy-five  percent 
greater  than  the  rates  predicted  with  the  one-dimensional 
model  with  linear  temperature  profile  along  the  wall.  This 
rather  significant  change  was  unexpected,  and  may  be  some¬ 
what  in  error  due  to  numerical  difficulties  in  the  method. 
These  numerical  difficulties  arose  due  to  the  boundary 
condition  at  node  1.  The  finite  element  boundary  conditions 
are  established  for  each  node  such  that  the  condition  is 
assumed  to  exist  over  the  entire  boundary  region  of  the  node. 
Therefore,  the  zero  flux  boundary  condition  could  not  be 
applied  over  the  whole  length  of  the  left  and  right  boundaries 
of  the  condenser  geometry.  This  problem  existed  for  nodes  1, 
4,  11,  and  12,  but  was  most  critical  for  node  1  since  it  was 
estimated  that  twenty  to  thirty  percent  of  the  total  heat 
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flux  through  the  inner  wall  passed  through  the  region 
associated  with  this  node. 

An  error  of  fifty  percent  in  the  calculation  of  node  1 
heat  transfer  rate  might  then  give  a  total  heat  transfer 
error  of  as  high  as  fifteen  percent.  Two  possible  solutions 
to  the  boundary  condition  problem  merit  attention.  The  first 
is  to  increase  the  number  of  elements  used  in  the  solution  of 
the  geometry  of  Figure  6.  However,  a  boundary  region  associ¬ 
ated  with  node  1  would  still  exist  on  the  left  hand  side  of 
the  geometry  although  reduced  in  size.  The  effects  of  the 
boundary  condition  problem  would,  therefore,  be  reduced  but 
not  eliminated. 

The  second  solution  would  be  to  consider  a  wall  geometry 
as  shown  in  Figure  15.  The  zero  flux  boundary  condition 
would  then  be  forced  along  the  boundary  in  question  due  to 
the  symmetry  of  the  problem.  The  use  of  one  of  these  alter¬ 
nate  solution  methods  might  contain  the  added  benefit  of 
eliminating  the  numerical  difficulties  associated  with  the 
stainless  steel  condenser  solution. 

The  heat  transfer  rate  solution  is  generally  affected  in 
the  same  manner  as  with  the  analysis  of  Schafer  when  h 

out 

0),  and  a  were  varied.  In  Figure  16  the  influence  of  outside 
convective  coefficient,  ,  is  shown.  Figure  17  illustrates 

how  the  rate  of  heat  transfer  was  relatively  unaffected  as 
fin  half  angle,  a,  was  varied.  In  addition,  heat  transfer 
rate  was  found  not  to  vary  significantly  with  rotation  rate 
as  shown  in  Figure  18.  This  result  concurred  with  the  results 
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obtained  earlier  by  Schafer  [8] .  Schafer  suggested  that  the 
invariance  in  heat  transfer  rate  with  the  rotation  rate  was 
due  to  the  cancelling  effects  of  decreasing  film  thickness 
and  decreasing  temperature  difference  across  the  film 
thickness  as  rotation  rate  is  increased. 

Because  of  numerical  difficulties  in  the  program,  no 
results  were  generated  for  the  stainless  steel  condenser. 

With  the  solution  of  the  first  and  second  test  cases,  the 
nodal  temperatures  calculated  were  accurate  only  to  within 
two.  or  three  degrees.  In  the  condenser  problem,  this  slight 
inaccuracy  resulted  in  nodal  temperatures  along  the  fin  wall 
which  were  slightly  above  the  fluid  saturation  temperature. 
Computer  output  indicated  that  the  wall  temperature  at  node  3 
was  approximately  0.1°F  above  the  saturation  temperature  when 
Tg  =  200°F.  As  a  result,  the  quantity  within  the  brackets  of 
equation  (23)  was  a  negative  value.  Error  messages  were 
generated  when  the  computer  attempted  to  evaluate  the  fourth 
root  of  this  negative  number. 
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IV.  CONCLUSIONS 


1.  The  addition  of  internal  axial  fins  to  the  condenser 
surface  of  a  rotating  heat  pipe  significantly  improves  heat 
transfer  performance, 

2.  Analysis  of  the  condenser  heat  transfer  problem  using 
a  two-dimensional  wall  conduction  model  and  a  parabolic 
temperature  distribution  along  the  fin  yielded  heat  transfer 
rates  about  seventy-five  percent  above  those  predicted  by  the 
one-dimensional,  linear  temperature  distribution  analysis  of 
Schafer  [8] . 

3.  For  the  geometrical  cases  studied  in  this  thesis,  the 
heat  transfer  rates  calculated  from  the  two-dimensional 
analysis  are  insensitive  to  changes  in  rotational  speed  and 
fin  half-angle. 

4.  Boundary  condition  problems  produced  with  the  finite 
element  solution  for  the  copper  condenser  were  the  probable 
cause  of  solution  inaccuracy  on  the  order  of  fifteen  percent. 

5.  Slight  inaccuracy  from  the  finite  element  calculation 
of  nodal  temperature  values  led  to  numerical  difficulties  in 
the  heat  transfer  solution  of  the  stainless  steel  condenser. 
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V.  RECOMMENDATIONS  FOR  FURTHER  STUDY 


It  is  suggested  that  the  following  topics  be  considered 
for  future  study  of  the  internally  finned,  rotating  heat 
pipe : 

1.  Re-examination  of  the  finite  element  analysis  using 
either  an  increased  number  of  elements  for  the  geometry 
considered  in  this  study  or  the  mirror  image  geometry  pre¬ 
sented  in  Figure  15. 

2.  Use  of  twelve-nodal  elements  in  the  finite  element 
analysis  to  allow  the  assumption  of  a  cubic  inner  wall  temper¬ 
ature  profile  in  order  to  more  closely  approximate  the  actual 
temperature  distribution. 

3.  Analysis  of  film  condensation  on  a  variety  of  fin 
shapes,  and  application  of  the  solution  to  the  two-dimensional 
model  to  determine  the  shape  permitting  optimum  heat  transfer 
rate. 

4.  An  experimental  program  in  order  to  compare  theoretical 
and  experimental  heat  transfer  rates  of  the  internally  finned 
condenser. 
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Figure  1.  Rotating  Heat  Pipe  Components 


37 


R 


Figure  2a.  Internally  Finned  Condenser  Geometry 
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I  r  =  R  +  X  sin  <}> 
1  2 
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Figure  2b.  Planar  Conduction  Simplification 
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Figure  3.  Condenser  Geometry  Considered  in  Analysis 
of  Schafer  [8] 
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R,  =  0.762  in . 
<D  =  1° 


Figure  4. 


Condenser  Geometry  Considered  in  Analysis 
of  Schafer  [8] 
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Figure  5.  Computer  Program  Flowchart  for  the  Analysis 
of  Schafer  [_8] 
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t  =  0,03125  in. 
b  =  0.025  in . 

R,  =  0.762  in . 

(f.  =  1° 


Figure  6,  Condenser  Geometry  Considered  in 
Two-Dimensional  Analysis 
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q  =  91822  Btu/Hr 
[91823] 


Figure  7.  Test  Case  (1) 
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T  ,  =  100°F 

oo  1 

hj  =  100  Btu/Hr-ft2-°F 
Actual  T  =  24.1°F 


=  1000  Btu/Hr-ft2-°F 
=  0°F 

q  =  7586  Btu/Hr 

[7527] 


Figu  re  8.  Test  Case  (2) 
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T  =  200°F 

h7  =  100  Btu/Hr-ft^-°F 


0. 75-c 

-r —  ft  - n 

hj  =  1000  Btu/Hr-ft^-°F 

T  ,=  100°F 
00  ^ 


C  N. 

One- 

Dimensional 

Two- 

Dimensional 

0.25  ft 

6526 

6727 

0.05  ft 

6561 

6754 

Figure  9. 


Test  Case  (3) 
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T  ,  =  200®F 

OO  1 

=  100  Btu/Hr-£t2 -°F 


,  0.2  5 _ _  0.5  _ _ 

r~  ft  ft  n 

=  1000  Btu/Hr-ft2-°F 
T  =  100°F 

00  Z 

q  =  6563  Btu 

[6921] 


Figure  10.  Test  Case  (4) 
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T  ,  =  200°F 
00  1 

hi  =  100  Btu/Hr-£t^-°F 


0.25  0.5 

ft  T  ft  n 

=  1000  Btu/Hr-ft^-®F 
T  =  100°F 

oo  Z 

q  =  6336  Btu/Hr 
[6548] 


Figure  11.  Test  Case  (5) 
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Figure  12.  Computer  Program  Flowchart  of 

Two-Dimensional  Conduction  Analysis 
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Figure  13.  Node  6  Boundary  Reg 


Condenser  Heat  Transfer  Rate,  BTU/Hr  (x  10'^) 


Saturation  Temperature,  °F 


Figure  14.  Condenser  Heat  Transfer  Rate  vs. 
Saturation  Temperature 
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Figure  15. 


Alternate  Condenser  Wall  Geometry 
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Condenser  Heat  Transfer  Rate,  BTU/Hr  (x  10"*) 


Figure  16.  Condenser  Heat  Transfer  Rate  vs. 
Saturation  Temperature 
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Condenser  Heat  Transfer  Rate,  BTU/Hr  (x  10"^) 


Figure  17. 


Condenser  Heat  Transfer  Rate  vs. 
Saturation  Temperature 


Condenser  Heat  Transfer  Rate,  BTU/Hr  (x  10"^) 


h  1000  Hr-ft2-  °F 
a  =  30° 


Figure  18.  Condenser  Heat  Transfer  Rate  vs. 
Rotation  Rate 
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FINC  TFE  NEXT  ITERATE  CF  TROUGH  CCNCENSATE  THICKNESS 
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